Possible d° ferromagnetism in MgO doped with nitrogen 
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We study the possibility of d° ferromagnetism in the compound MgO doped with nitrogen (N). 
The Haldane- Anderson impurity model is formulated within the tight-binding approximation for 
determining the host band-structure and the impurity-host hybridization. Using the quantum Monte 
Carlo technique, we observe a finite local moment for an N impurity, and long-range ferromagnetic 
correlations between two N impurities. The ferromagnetic correlations are strongly influenced by the 
impurity bound state. When the ferromagnetic correlation between a pair of impurities is mapped 
onto the isotropic Heisenberg model for two spin- 1/2 particles, the effective exchange constant 
J12 is found to increase with increasing temperature. Similar temperature dependence of J\% is 
also obtained in other diluted magnetic semiconductors, such as zincblende ZnO doped with Mn. 
The temperature dependence of J12 suggests that the mapping of the full Hamiltonian onto the spin 
Hamiltonian cannot fully describe the magnetic correlations for the diluted magnetic semiconductors 
at least in the limit of low impurity spin. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.40.Mg 



I. INTRODUCTION 



There are currently a large number of observations 
of ferromagnetism associated with doping of a pri- 
ori non-magnetic species into non-magnetic oxide semi- 
conductors [l|, H, The subject has been given the 
name of d° magnetism to emphasize the fact that the 
magnetism is probably not coming from partially filled 
d-orbitals, but from moments induced in the p-orbitals 
of the oxygen band 0, [H, @|. Even in cases where 
there are partially filled d-orbitals, discrepancies between 
bulk measurements indicating ferromagnetism, and mi- 
croscopic measurements by XMCD finding paramagnetic 
transition metal ions, have suggested that the ferromag- 
netism comes from vacancies in the oxygen lattices [7|. 
While a fertile area of current investigation, the subject 
remains obscure because of problems of irreproducibil- 
ity for different samples, sample inhomogeneity, stability 
in time, etc., and the exact origin of the magnetism is 
not yet established. This is unlike the case of regular di- 
luted magnetic semi-conductors such as (Ga,Mn)As Q, 
where ferromagnetism is reproducible in controlled sam- 
ples. The theoretical situation is also much less clear for 
d° magnetism than in the case of "classic" diluted mag- 
netic semiconductors where the origin of the magnetic 
moments, the substituted magnetic ion, is in no doubt 
and, in the case of Mn 2+ , quite large (S=5/2). Mag- 
netic exchange interactions leading to ferromagnetism 
are mediated by holes in the p band. For (Ga,Mn)As, 
Curie temperatures can be estimated with reasonable ac- 
curacy, starting from band structure to first derive the 
effective exchange at different distances between spins in 
an effective spin Hamiltonian 0] which is then studied 
via classical Monte Carlo or other approximate methods 
pLllllli3|. Although this approach is adequate for bulk 
properties such as the ordering temperature for a variety 
of materials, fuller treatment of the effects of the corre- 



lation are needed for understanding more detailed prop- 
erties such as the local densities of states and impurity 
band structure (l3| . Even for thermodynamic proper- 
ties such as the Curie temperature, there is good reason 
to doubt straightforward application of methods appar- 
ently working for "classic" diluted magnetic semiconduc- 
tors to the materials that may show d° magnetism. In 
this case, the induced moments postulated in the oxy- 
gen bands are not fully localized, but are formed by the 
same holes that generate effective exchange. It was re- 
cently argued, for example, that [l4[ because of this, 
exchanges estimated from local approximations, such as 
the local spin density approximation (LSDA) may con- 
sistently over-estimate the tendency for ferromagnetism 
and even formation of local moments, compared to ex- 
tensions, such as self-interaction correction (SIC). Thus 
much of the theoretical literature may be over-optimistic 
in predictions of ferromagnetism, even at zero temper- 
ature. One can even question whether the approach is 
valid to first project the full Hamiltonian onto the zero- 
temperature spin Hamiltonian, to which the effects of 
thermal fluctuations are added. 

It is important, therefore, to develop methods which 
treat correlations correctly and that do not rely on arbi- 
trary approximations or on the separation between spin 
and charge fluctuations. Such a method has been in- 
troduced and applied to models of the "classic" diluted 
magnetic semiconductors [TEl . [l6l . Il7| , namely Quantum 
Monte Carlo methods based on the Hirsch-Fye algorithm. 
These can directly give spin correlations at finite tem- 
peratures without making any assumption of projection 
onto an effective Hamiltonian. In this paper we apply the 
method to an interesting case of possible d° magnetism, 
namely MgO diluted with nitrogen. Our aim is then to 
answer, by means of Quantum Monte Carlo methods the 
following questions. Firstly, does such an unbiased calcu- 
lation predict ferromagnetic correlations in such a mate- 
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rial? Secondly, of more general interest, do the standard 
approaches, as described above, still apply? In fact fer- 
romagnetism has previously been predicted at low con- 
centrations [IH in the doped alkaline-earth metal oxides 
(MgO, CaO, BaO, SrO) doped with N and C. Our re- 
sults are not completely comparable to those results, in 
that we shall consider a simpler model Hamiltonian for 
the host and two impurities only, rather than a finite 
concentration. We should be able to draw useful gen- 
eral conclusions both for the particular material and the 
methods. 



II. FERROMAGNETIC CORRELATIONS 
BETWEEN TWO N IMPURITIES 

A. Impurity model 

In order to describe N impurities in MgO host, we take 
the two-step calculations. Firstly, the Haldane- Anderson 
impurity model [T|| is formulated within the tight- 
binding approximation for determining the host band 
structure and the impurity-host hybridization. Secondly, 
the magnetic correlations of the Haldane-Anderson im- 
purity model at finite temperatures are calculated by the 
Hirsch-Fye quantum Monte Carlo technique [20| . 

The Haldane-Anderson impurity model is defined as 

k,a,cr 

+ X! (^i?kQPi ?cr Ck QC r + H.C.) 
k,a,i,£ ,ct 

where c^ a(7 (ck Q(T ) is the creation (annihilation) operator 
for a host electron with wavevector k and spin a in the 
valence (a = v) or conduction (a = c) band, and p\^ a 
(Pi(a-) is the creation (annihilation) operator for a local- 
ized electron at impurity site i in orbital £ (£ = x, y, 
z) and spin a with m^ a — p\^ a Pie,<j- Here, e Q (k) is the 
host band dispersion, fx the chemical potential, V^ka the 
mixing between the impurity and host, e p the impurity 
2p orbital energy, and U the on-site Coulomb repulsion 
for the impurity. 

The energy bands e Q (k) for MgO host, and the 
impurity-host hybridization Vi^ka will be calculated 
within the tight-binding approximation. For a large num- 
ber of simple oxides, the on-site Coulomb repulsion en- 
ergy of holes in an oxygen p orbital is 5-7 eV [21| . For the 
on-site Coulomb repulsion of 2p-orbitals at an N impu- 
rity site in the compound MgO, the experimental value 
is unknown, so here we take it as U = 6 eV. In addi- 
tion, the experimental value of impurity 2p energy e p (N) 
in Mg(0,N) is also unknown. In the following, we use 
e p (N) as the parameter satisfying e p (N) > e p (0)= -2.03 
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FIG. 1: (Color online) Energy bands of MgO with rocksalt 
crystal structure. These results were reproduced using the 
tight-binding parameters in Ref. [S^], where one 3s orbital of 
Mg and three 2p orbitals of O are included. 



eV, where the value of e P (0) is taken from the tight- 
binding parameters for the compound MgO (see below). 

For the compound Mg(0,N), there is one 2p hole at an 
N 2 ~ site, so it is reasonable to neglect the Hund couplings 
among different 2p orbitals at an N 2_ site. 



B. Tight-binding approach for the MgO band 
structure and the N-MgO hybridization 

In this section, we discuss the tight-binding calcula- 
tion for the MgO band-structure, and the hybridization 
between an N impurity and MgO host. For MgO with 
rocksalt structure, the band structure e Q (k) had already 
been calculated with a set of tight-binding parameters 
[22l |. In this approach, the basis consists of a 3s orbital 
for the cation Mg 2+ and three degenerate 2p orbitals for 
the anion 2 ~. The orbital energies are e s (Mg) = 9.88 
eV, e p (0)= -2.03 eV. In addition, the mixing values be- 
tween the s orbital of Mg 2+ and the p orbitals of 2 ~ 
up to the 3rd- nearest neighbors are taken to be (sper)i 
= 1.10 eV, (ssa) 2 = -0.18 eV, {ppa) 2 = 0.65 eV, (pptt) 2 
= -0.07 eV, and (spa) 3 = 0.89 eV, where () n means the 
mixing of the n-th nearest neighbors. Using these tight- 
binding parameters, we have reproduced the band struc- 
ture of rocksalt MgO as shown in Fig. [TJ where band 
structure consists of one conduction band and three va- 
lence bands. The conduction band mainly comes from 3s 
orbital of Mg, and valence bands mainly come from 2p 
orbitals of O. The top of valence bands and the bottom of 
conduction band are located at the T point with a direct 
gap of 7.72 eV. 

Next, we discuss the calculation of the hybridization 
between N impurity and MgO host within the tight- 
binding approximation. The hybridization matrix ele- 
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ment Vi^ka = (<Pz(i)\H l^a(k)) has the form of 

" o,n 

l 



where y^(i) is the impurity 2p-state (£ = x, y, z) at site i, 
and fy a (k) is the host state with wavevector k and band 
index a, which is expanded by atomic orbitals (/? ( n ) with 
orbital index o and site index n. Here, N is the total 
number of host lattice sites, and a QO (k) is an expansion 
coefficient. For the mixing integrals of ((p^(i)\H\(p (n)) , 
£ denotes the three 2p orbitals of N 2_ , and o represents 
the 3s orbital of Mg 2+ and 2p orbitals of 2 ~. As shown 
by Slater and Koster [23|], these mixing integrals up to 
the 3rd-nearest neighbors can be expressed by four inte- 
grals (spa)i (ppa)2, (pf>7r)2> (spa) 3 and direction cosines 
I, m and n in the two-center approximation, where ()„ 
means the mixing of the n-th nearest neighbors. As the 
experimental values of the above four integrals are un- 
known, here we take these mixing integrals between N 
and MgO as the same mixing values between O and MgO 
in Ref . [22| . Thus, we have {spa)\ = 1.10 eV, (j>po)i — 
0.65 eV, (ppir) 2 = -0.07 eV, and (spa) 3 = 0.89 eV for the 
impurity-host hybridization. 

Figure [2] displays results on the impurity-host mixing 
function V^(k) defined by 



V c (k)^(£|V^(k)| ; 



»V2 



(3) 



where £ is a 2p(x, y, z) orbital of an N impurity. In 
Eq. ([3]), the summation over a is performed only over 
the valence bands (Fig. EJa)) or the conduction band 
(Fig. IHb)). Here, V^(k) is plotted along various direc- 
tions in the Brillouin zone for rocksalt crystal structure. 
It is observed that, at the T point, the total hybridiza- 
tion between the £ orbital of an N impurity and the MgO 
valence bands is finite, while that with the MgO conduc- 
tion band is zero. For the host MgO, the gap edge is 
located at the T point, hence the value of V near T will 
be particularly important in determining the energy of 
the impurity bound state, which may appear in the gap 
due to the mixing between impurity and host. 



C. Quantum Monte Carlo results 

In this section, we present results on the magnetic cor- 
relations of the Haldane-Anderson impurity model ob- 
tained by the Hirsch-Fye quantum Monte Carlo (QMC) 
technique [2(| • The parameters related to the MgO host 
and N impurity have been calculated within the tight- 
binding approach described above. The following results 
were obtained with more than 10 5 Monte Carlo sweeps 
and Matsubara time step At = 0.225. 
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FIG. 2: (Color online) Hybridization function V^(k) between 
£ orbital of an N impurity and MgO host (a) valence bands or 
(b) conduction band. Here, the mixing parameters between 
N and MgO are taken as the same mixing values between O 
and MgO in Ref. [13]. 



We first discuss the local moment formation for an N 
impurity in the MgO host. For this purpose, we have 
performed QMC simulations to calculate ((M z ) 2 }, where 



AI Z = n i4T - n ia 



(4) 



is the magnetization operator for a £ orbital of an N 
impurity at site i. Fig. [3] shows ((M z ) 2 ) versus the 
chemical potential /i at temperature T = 200K, where 
< fi < 0.5eV. As mentioned before, the experimen- 
tal value of impurity 2p energy e p (N) in Mg(0,N) is not 
known, so it is taken as a parameter satisfying the rela- 
tion e p (N) > e p (0) = -2.03 eV. As shown in Fig. H the 
sharp increases in the magnitude of the magnetic moment 
are observed around 0.08eV and 0.35eV for the impurity 
2p energy e p = — 1.5eV and e p = — 0.5eV, respectively. In 
our model 2p orbitals are degenerate, and thus the calcu- 
lated curves in Fig. [3] do not change with £ = x, y or z. 
In addition, no such sharp increases (or sharp decreases) 
are observed near the bottom of the conduction band. As 
shown in Fig. [21 the values of the hybridization with bot- 
tom of conduction band are around zero and thus much 
weaker than those with top of valence bands. 

According to the Hartee-Fock [24| and QMC calcula- 
tions (TBI, [Ty, the presence of a sharp increase (or 
a sharp decrease) in ((M z ) 2 ) versus /i implies the exis- 
tence of an impurity bound state (IBS) at this energy, 
and the IBS plays an important role in determining the 
strength of the ferromagnetic (FM) correlations. When 
the IBS is unoccupied, FM correlations can develop be- 
tween the impurities. When the IBS is occupied, the FM 
correlations become weaker. 

When two N impurities are introduced in the MgO 
host, they can replace any two O positions in the rock- 
salt lattice MgO. Let us then briefly explain how our 
QMC calculations of the magnetic correlations function 
(Mf M|) actually proceed with the different spatial po- 
sitions of the two impurities. When the position of the 
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FIG. 3: (Color online) Square of magnetic moment ((M z ) 2 ) 
of an N impurity versus chemical potential [i at T = 200K. 
The experimental value of impurity 2p energy e P (N) in 
Mg(0,N) is unknown, but it should satisfy the relation e p (iV) 
> tp{0) = —2.03 eV. Here, the arrows indicate the different 
position of fi investigated by Fig. [4] 



two impurities changes, the corresponding hybridization 
Vl^ka ; including the impurity spatial position i as defined 
in Eq.(2), will change; the Green's function between im- 
purity 1 and impurity 2, including the factor V*^ a V2^k a , 
will change; and thus the magnetic correlations function 
(MfMf) changes. See Ref.[25| for more details of the 
QMC calculations. 

The magnetic correlation function (MfMf) between 
£ — x orbitals of the impurities versus the impurity sep- 
aration R/a at temperature T = 200K is shown in Fig. 
IH The direction i?//(0.5, 0, 0.5) is chosen to be along 
one of the 12 nearest N-N neighbors in rocksalt struc- 
ture, and a is the lattice constant. Results in Fig. [4ja) 
are obtained with fixed impurity energy e p = —1.5 eV for 
various chemical potential \x. The impurity spins exhibit 
FM correlations at chemical potential fi = 0.0 eV. It is 
noted that the IBS of Mg(0,N) lies near 0.08 eV in the 
gap, close to the top of the valence band as displayed in 
Fig. [3J By increasing to 0.06 eV, the FM correlations 
become larger, and the range of the FM correlations be- 
comes longer. Further increasing /i, the FM correlations 
become weaker. This is because the IBS of Mg(0,N) be- 
comes occupied when \i is increased to above 0.06 eV, as 
seen in Fig. [3] Results in Fig. [4jb) are obtained with 
fixed impurity energy e p = —0.5 eV for various values 
of the chemical potential \x. The impurity spins exhibit 
quite weak antiferromagnetic (AFM) correlations due to 
the superexchange interaction at fi = 0.0 eV. It is worth 
pointing out that the IBS of Mg(0,N) lies around 0.35 
eV, deep in the band gap as shown in Fig. [3] Increasing 
[i to 0.3 eV, the FM correlations appear with stronger 
magnitude and longer range. Further increasing fi to 
above 0.3 eV, the FM correlations become weaker. This 
is because the IBS of Mg(0,N) is occupied when /i > 0.3 
eV, as displayed in Fig. [3] 
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FIG. 4: (Color online) Magnetic correlation function 
{Mi Mi) between £ = x orbitals of the N impurities versus 
distance R/a for (a) the impurity 2p energy e p = -1.5 eV and 
(b) e p = -0.5 eV at T = 200A'. Direction R//(0.5, 0, 0.5) is 
along one of the 12 nearest N-N neighbors in rocksalt struc- 
ture. 



In the direction i?//(0.5, 0, 0.5), let us discuss the cor- 
relation function (Mf M|) between other orbitals of two 
N impurities in the MgO host. Considering the sym- 
metry in the direction R//(0.5, 0, 0.5), (Mf Mf ) between 
£ = z orbitals of two N impurities have the same value 
between £ = x orbitals of two N impurities. For the 
(MfMf) between £ = y orbitals of two N impurities in 
the direction R/ / (0.5, 0, 0.5), our QMC calculations (not 
present here) show that it is short-range FM correlation. 
For the nearest N-N neighbors distance R = (0.5, 0, 0.5)a 
in the case that e p = -1.5 eV, fi = 0.06 eV and T = 200K, 
the magnetic correlation function (MfMJ) between dif- 
ferent £ orbitals is shown in Fig. [5l It is found that the 
magnetic correlations between different orbitals of two N 
impurities are much smaller than the values between the 
same orbitals. 

For the two N impurities in the MgO host, there 
are 12 nearest N-N neighbors in rocksalt structure : 
(0,±0.5,±0.5)a, (±0.5,0,±0.5)a, (±0.5, ±0.5, 0)a. For 
the direction R along other 11 nearest N-N neighbors, 
(MfM|) of two N impurities could be got by sym- 
metry without calculation. For example, (MfM|) be- 
tween £ = x orbital of two N impurities in direction 
R//(0, 0.5, 0.5) is the same value as that between £ = y 
orbitals of two N impurities in direction R//(0.5, 0, 0.5). 
Thus, for e p = -1.5 eV and fi = 0.06 eV, the long-range 
FM correlation between two N impurities could be ob- 
served at T = 2007^ along all of the 12 nearest N-N 
neighbors in rocksalt Mg(0,N). 
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FIG. 5: (Color online) Magnetic correlation function 
(MfMf) between 2p orbitals of two N impurities. Distance 
R — (0.5, 0, 0.5)a is one of the 12 nearest N-N neighbors in 
rocksalt structure. Here, e p — -1.5 eV, \i = 0.06 eV and 
T = 20074". 



The calculations presented above do not include the 
Jahn-Teller(JT) distortion, which might occur in the 
compound Mg(0,N). As the JT distortion is included, 
the energy level of x, y, z orbitals of the impurity N will 
differ, being higher for some orbitals and being lower for 
the others, but the values changed are often quite small. 
As a result, the IBS shown in Fig. [3] will differ for differ- 
ent p orbitals, being shallower in the gap for those having 
lower energy level and being deeper in the gap for those 
having higher energy level, and again the shift values of 
the IBS will be quite small. More importantly, the long- 
range FM correlations, as shown in FigHJ will become 
stronger for those having the shallower IBS and become 
weaker for those having the deeper IBS, and of course the 
change will be quite small. So, here we argue that the JT 
distortion will play a quite small role in our calculations. 



III. MAPPING OF FERROMAGNETISM ONTO 
A HEISENBERG MODEL 

The ferromagnetic correlation between two N impuri- 
ties is mapped onto the isotropic Heisenberg model for 
two spin- 1/2 particles 



H = - Ji 2 Mi • M 



(5) 



At finite temperature, defining j3 = l/k B T with k B the 
Boltzmann constant, the impurity-impurity correlation 
is defined as 



(MfMf) 
Considering 

Mi • M 2 



Tr(M?M£e-P H )/Tr(e-P H ). 



;(M! + m 2 ) 2 



M 2 



Mi 



where s = 1/2, S = 0,1, and the trace here could be 
taken as 



Tr< 



we have 



E £(• 

s z =~s s=o 



1/2 1/2 

E E (•••), 

s |=-l/2sf=-l/2 



(Mi MS ) = -■ - — € -^— ■ 

\ 1 2 / 4 3 + e -/3Ji2' 

where the unit of (MfMf) is l 2 = (2s) 2 = (2/i B ) 2 , fi B 
is Bohr magneton. To be consistent with our QMC cal- 
culation results, whose unit of (MfMf) is /i B , the above 
equation is modified as 



(MfM| ) = 



3 + e~P J ™ ' 



(6) 



Thus, inverting this relation, we can deduce an effective 
exchange coupling J12 between two N impurities from 



J12 = k B T In 



1 + (Mf mi; 

1 - 3 (MfMf 



(7) 



S(S+l)-2s(s+l)], 



Figure [6J shows the temperature dependence of the 
magnetic correlation function (MfMf) and the corre- 
sponding exchange coupling J12 between £ = x orbitals 
of two N impurities. Here, the following parameters are 
taken as e p = -1.5 eV, fi = 0.06 eV, R = (0.5,0,0.5)a 
or R — (1,0, l)a. For the two N impurities with dis- 
tance R = (0.5, 0, 0.5)a, it is found that with increasing 
temperature from 20QK to 400K, magnetic correlation 
function (MfMf) decreases, but the exchange coupling 
Ji 2 increases. For the two N impurities with distance 
R = (1,0, l)a, the temperature dependence of exchange 
coupling J12 becomes weaker. 

For (MfMf) and J12 between £ = y or z orbitals of two 
N impurities in the direction i?//(0.5, 0, 0.5), the similar 
behaviors are observed. For the direction R along other 
11 nearest N-N neighbors, the behaviors of temperature 
dependence could be found by symmetry without calcu- 
lation as discussed in the last section. Thus, with e p = 
-1.5 eV and \x = 0.06 eV, the temperature dependence of 
(MfMf) and J12 between two N impurities, as shown in 
Fig. [51 could be observed along all of the 12 nearest N-N 
neighbors in rocksalt Mg(0,N). 



IV. DISCUSSION 

To understand the long-range FM correlation func- 
tion (MfMf) between two N impurities in MgO host, 
as shown in Fig. [H which is mediated by the impurity- 
induced polarization of the host electron spins, it is more 
convenient to study the local density of states in real 
space around the impurity. Here we have considered the 
impurity-host correlation function (M z m z (r)) and the 
number of polarized host electrons (n(r) — n(oo)) for the 
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FIG. 6: (Color online) Temperature dependence of (a)the 
magnetic correlation function (Mf Af| ) between x orbitals of 
two N impurities, and (b) corresponding exchange coupling 
J12 calculated by Eq.©. Here, e p = -1.5 eV, ^ = 0.06 eV, 
R = (0.5, 0, 0.5)a or R = (1, 0, l)a. 
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FIG. 7: (Color online) Spatial variation of (a) the magnetic 
correlation function {M z m z (r)) between x orbital of an N 
impurity and MgO host, and (b) the number of polarized 
host electrons (n(r) — n(oo)). Here, e p = -1.5 eV, \i = 0.06 
eV, and r//(0.5, 0,0.5). 



case of one- impurity N in the MgO host. Here, r is the 
site of host electron and the impurity N is located at site 
r = 0. In addition, n(oo) means the number of host elec- 
trons at infinity distance, thus (n(r) — n(oo)) is negative 
and represents the number of holes. The magnetization 
m z (r) and number n(r) operators are defined as 

™ Z ( r ) = ^2{n a v] - n arl ), (8) 

a 

n ( r ) = ^2(n a ri +n arl ), (9) 

a 

where n arcr — c? aTu c a ra is the number operator for host 
electrons with band index a and site r and spin a. With 
impurity 2p energy e p = -1.5 eV and chemical potential 
(i = 0.06 eV and direction r//(0.5, 0, 0.5), the long-range 
antifcrromagnetic (AFM) correlation between x orbital 
of N impurity and MgO host is observed in Fig. [Tfa), 
and the polarized host electrons with long-range distri- 
bution are also observed in Fig. \T[b) . As the holes in the 
valence band round the impurity N are spin-polarized, 
the valence band is spin-polarized. Comparing Fig. [3^a) 
with Fig. [Zja), it is confirmed that the long-range 
AFM impurity-host correlations contribute to the long- 
range FM impurity-impurity correlations. With increas- 
ing temperature, as shown in Fig. [7(a), the magnitude 
of AFM impurity-host correlation function (M z m z (r)) 
decreases, which induces the decreasing FM impurity- 
impurity correlation function (MfM|) with increasing 
temperature as shown in Fig. [5(a). 



To understand the temperature dependent exchange 
coupling J i2 shown in Fig. EJb), we have studied 
a "classic" diluted magnetic semiconductor, zincblende 
ZnO doped with Mn. As shown in Fig. [51 a simi- 
lar temperature-dependent exchange coupling J12 is ob- 
tained. Here, magnetic correlation function (MfMf) 
between xy orbitals of two Mn impurities is calculated, 
and the corresponding exchange coupling J12 is given by 
Eq.([7j). Because the experimental value of for Mn in 
ZnO host is unknown, here we use the symmetric case 
of €d — -U/2 + /i so that the impurity sites develop 
large magnetic moments. For the compound (Zn,Mn)0, 
the value of the on-site Coulomb repulsion for Mn 2+ is 
taken as U — 5.2eV by comparing with the photoemission 
spectroscopy measurements [26| . The chemical potential 
value fj, = 0.1 eV is set close to the impurity bound state. 
The direction RJ / (0.5, 0, 0.5) is chosen to be along one of 
the 12 nearest Mn-Mn neighbors in zincblende structure. 
We find that the exchange constant J12 increases with 
increasing temperature as was the case in Mg(0,N). 

We have to note that our result of J12 does indeed con- 
tradict the common assumption that people make when 
using the Lichtenstein formula, but here we are dealing 
with a low spin situation, in contrast to what is usually 
(but not always) done in the literature. Thus our results 
are most pertinent for d° situations, and the results for 
zincblende ZnO doped with Mn is in fact based on an 
effective spin-1/2 model. 



7 



A 

N <M 



N i- 

V 



0.1 



(a) 



(Zn,Mn)0, ^ = xy, |i = 0.1 eV 



R = (0.5, 0, 0.5)a — b- 
(1,0, 1)a 








i 




300 


_(b) 


i 




CD 
J£ 


200 










100 













i 



200 300 400 

Temperature (K) 

FIG. 8: (Color online) For zincblende ZnO doped with Mn, 
temperature dependence of (a) the magnetic correlation func- 
tion (MfM|) between xy orbitals of two Mn impurities, and 
(b) the corresponding effective exchange coupling J12 calcu- 
lated by Eq. (0. Here e d = -U/2 + fj,, (1 = 0.1 eV, and 
R//(0.5, 0,0.5). 



V. SUMMARY AND CONCLUSIONS 

In summary, we have studied possible d° ferromag- 
netism for the compound Mg(0,N) in the dilute impurity 
limit based on the Haldane- Anderson impurity model. 
The band structure of the MgO host were calculated us- 
ing the tight-binding approach. The mixing parameters 



between N and MgO are approximated as being the same 
as between O and MgO. The QMC results show the de- 
velopment of a large magnetic moment at an N impurity 
site, and long-range ferromagnetic correlations between 
two N impurities. The ferromagnetic correlation between 
impurity pairs is mapped onto the isotropic Heisenberg 
model for two spin- 1/2 particles, and the effective ex- 
change coupling J12 for given separation of impurities is 
found to increase with increasing temperature. Similar 
temperature dependence of J12 is also obtained in "clas- 
sic" diluted magnetic semiconductors, such as zincblende 
ZnO doped with Mn, suggesting that the mapping is not 
fully valid even there at least in the limit of low impurity 
spin. For the particular case of MgO doped with N the 
results presented in this paper, in which no such approx- 
imation is made and in which interactions are treated 
exactly, suggest that there should be stable moments as- 
sociated with each impurity and that they should have 
ferromagnetic correlations. While the results are for two 
impurities only, the long range of the correlations sug- 
gests long-range ferromagnetic order for low concentra- 
tions. 

ACKNOWLEDGMENTS 



This work was supported by the NAREGI Nanoscience 
Project and a Grant-in Aid for Scientific Research from 
the Ministry of Education, Culture, Sports, Science and 
Technology of Japan, and NEDO. The authors thank the 
Supercomputer Center at the Institute for Solid State 
Physics, University of Tokyo, for the use of the facilities. 
T. Z. thanks the International Frontier Center for Ad- 
vanced Materials of Tohoku University and the members 
of the IMR for their support of his stay in Sendai, which 
enabled this collaboration. The authors acknowledge S. 
Parkin for valuable discussion about the experiment of 
Mg(0,N). 



[1] M. Venkatesan, C. B. Fitzgerald, and J. M.D. Coey, Na- 
ture (London) 430, 630 (2004). 

[2] M. Venkatesan, C. B. Fitzgerald, J. G. Lunney, and J. 
M. D. Coey, Phys. Rev. Lett., 93, 177206 (2004). 

[3] J. M. D. Coey, M. Venkatesan, P. Stamenov, C. B. 
Fitzgerald, and L. S. Dorneles, Phys. Rev. B72, 024450 
(2005). 

[4] A.M. Stoneham, A.P. Pathak and R.H. Bartram, J. Phys. 

C: Solid State Phys. 9, 73 (1976). 
[5] I. S. Elfimov, S. Yunoki, and G. A. Sawatzky, Phys. Rev. 

Lett. 89, 216403 (2002). 
[6] C. Das Pemmaraju and S. Sanvito, Phys. Rev. Lett. 94, 

217205 (2005). 

[7] T. Tietze, M. Gacic, G. Schtz, G. Jakob, S. Brck and E. 

Goering, New J. Phys. 10 055009 (2008). 
[8] H. Ohno, Science 281,951 (1998). 

[9] A. I. Liechtenstein, M. I. Katsnelson, V. P. Antropov, and 
V. A. Gubanov, J. Magn. Magn. Mater. 67, 65 (1987). 
[10] L. Bergqvist, O. Eriksson, J. Kudrnovsky, V. Drchal, P. 
Korzhavyi, and I. Turek, Phys. Rev. Lett. 93, 137202 



(2004). 

[11] K. Sato, W. Schweika, P. H. Dederichs, and H. 

Katayama-Yoshida, Phys. Rev. B 70, 201202(R) (2004). 
[12] G. Bouzerar, T. Ziman and J. Kudrnovsky, Europhys. 

Lett., 69, 812-818 (2005). 
[13] J. Ohe, Y. Tomoda, N. Bulut, R. Arita, K. Nakamura, 

and S. Maekawa, unpublished (arXiv.0812.0430). 
[14] A. Droghetti, C. D. Pemmaraju and S. Sanvito, Phys. 

Rev. B 78, 140404(R) (2008). 
[15] N. Bulut, K. Tanikawa, S. Takahashi, and S. Maekawa 

Phys. Rev. B 76, 045220 (2007). 
[16] Y. Tomoda, N, Bulut, and S. Maekawa, to be publised 

in Physica B (!arXiv:0806.0095|l . 
[17] B. Gu, N. Bulut, and S. Maekawa, J. Appl. Phys. 104, 

103906 (2008). 

[18] K. Kenmochi , V. A. Dinh, K. Sato, A. Yanase, H. 

Katayama-Yoshida, J. Phys. Soc. Jpn.73, 2952 (2004); V. 

A. Dinh, M. Toyoda, K. Sato, and H. Katayama-Yoshida, 

J. Phys. Soc. Jpn.75 093705 (2006). 
[19] F.D.M. Haldane and P.W. Anderson, Phys. Rev. B 13, 



8 



2553 (1976). 

[20] J. E. Hirsch and R. M. Fye, Phys. Rev. Lett. 56, 2521 
(1986). 

[21] I.S. Elfimov, S. Yunoki and G.A. Sawatzky, Phys. Rev. 

Lett. 89, 216403 (2002). 
[22] V. C. Lee, and H. S. Wong, J. Phys. Soc. Jpn. 45, 895 

(1978). 

[23] J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 
[24] M. Ichimura, K. Tanikawa, S. Takahashi, G. Baskaran, 



and S. Maekawa, Foundations of Quantum Mechanics 

in the Light of New Technology, eds. S. Ishioka and K. 

Fujikawa, (World Scientific, Singapore, 2006), 183-186, 

( |cond-mat/0701736j ). 
[25] R. M. Fye, J. E. Hirsch, and D. J. Scalapino, Phys. Rev. 

B 35, 4901 (1987). 
[26] T. Mizokawa, T. Nambu, A. Fujimori, T. Fukumura, and 

M. Kawasaki, Phys. Rev. B 65, 085209 (2002). 



